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Abstract. We develop a formalism for the calculation of the flow of angular 
momentum carried by the fluctuating electromagnetic field within a cavity 
bounded by two flat anisotropic materials. By generalizing a procedure employed 
recently for the calculation of the Casimir force between arbitrary materials, we 
obtain an expression for the torque between anisotropic plates in terms of their 
reflection amplitude matrices. We evaluate the torque in ID for ideal and realistic 
model materials. 
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1. Introduction 

In the last decade the Casimir effect pQ has received considerable attention, as 
the recently attained high experimental accuracy has permitted detailed tests of 
theoretical predictions [21|310I|Slini|71|Slin|- This, in turn, has stimulated a 
growing interest in fundamental aspects of the vacuum field. The study of vacuum 
forces between realistic materials was pioneered by Lifshitz ^U], who considered local 
homogeneous materials whose fluctuating currents were the sources of the fluctuating 
electromagnetic field and whose correlations were related to the dielectric response 
of the materials. One of the limitations of the Lifshitz theory is the requirement of 
a definite microscopic model of matter which has to be solved simultaneously with 
the electromagnetic field equations. Thus, the applicability of the results seem to 
be limited by the generality of that initial model. In particular, Lifshitz results were 
developed for isotropic materials. Nevertheless, the vacuum energy has been calculated 
for cavities bounded by anisotropic materials, first in the non-retarded limit and 
later for arbitrary distances jj^lj resulting in a torque whenever the optical axes of 
the plates are not aligned with each other. An alternative derivation of the Casimir 
torque in the ID case has been developed starting from the angular momentum 
flux carried by the fleld. Analytical formulae have been found in the retarded limit 
when the anisotropy is small |14| . Numerical calculations have also been performed for 
materials with a small anisotropy and it has been shown that the torque may be large 




Figure 1. (a) Vacuum cavity V of width L bounded by two arbitrary material 
slabs (1 and 2) with surfaces at zi and Z2. A photon with polarization 
propagating within the cavity is either reflected coherently with amplitude and 
polarization r2e" = r^^el^, or lost from the cavity with probability 1 — |rp. 
(b)Fictitious system made up three empty regions /, II, and bounded by 

perfect mirrors at zq and 23 and with infinitely thin sheets at 21 and 22 whose 
reflection amplitudes r"'' are identical to those of the real system and whose 
transmission amplitudes t"'' are such that energy is conserved with no absorption 
whatsoever. The field within the fictitious cavity V' (region //) is the same as 
within the real cavity V. 



enough to be experimentally measurable in several novel experimental configurations 

In the previous theoretical works essential assumptions about the dielectric 
properties of the plates were done from the onset in order to derive expressions for 
the Casimir torque. However, recent works \W\ I17| have shown that if the theory is 
set up in terms of the reflection coefficients of the media, it is possible to decouple 
the calculation of the Casimir force from the calculation of the dielectric response 
of the materials. The so called scattering approach has permitted the calculation of 
the force for a wide class of systems simply by plugging into the resulting formulae 
the appropriate reflection amplitudes or surface impedances. Thus, transparent and 
opaque, local and non-local, infinite and finite, homogeneous and heterogeneous 
systems may be treated in the same footing. In the present paper we generalize the 
scattering approach to account for anisotropy as well. We present a new derivation of 
the Casimir torque between plates with arbitrary dielectric properties characterized 
by their anisotropic optical coefficients. For simplicity, in this paper we focus our 
attention on one dimensional systems, although our approach is also applicable to 
3D PHI- We present results for both ideal systems, for which analytical formulae are 
obtained, and realistic dichroic systems. 

2. Scattering approach 

To calculate the torque we follow the scattering approach I17| . illustrated by 
Fig.^a). A photon within the cavity may be described by its amplitude f , frequency 
uj, wavevector k = (Q, ±g) and polarization e". We may chose the independent 
polarizations as a = s,p (or equivalently, TE and TM). When a photon with 
polarization ef is incident upon the surface, say, of medium 2 at Z2, it is either 
reflected coherently with an amplitude and polarization proportional to r2e" = rj'^ef 
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(sum implied), or transmitted with a probability 1 — |7'2p- Here, r"^ is a 2 x 2 matrix 
that describes the reflection amplitude of medium a = 1,2. For isotropic media, r"^ 
may be taken as a scalar whenever the incoming field has s or p polarization, but 
that separation is not possible in general when the media are anisotropic. We remark 
that r^^ is defined to be the complete reflection amplitude of medium a, not only 
that of its front surface. For instance, if medium a were a thin film or a layered 
system, the multiple reflections within a are to be incorporated into r^^ . Thus, if 
the photon is not coherently reflected, it must necessarily be absorbed within a or 
else be transmitted into the empty space beyond. The principle of detailed balance 
implies that in thermodynamic equilibrium, for every photon that is not reflected 
and is therefore lost from the cavity either through absorption or transmission, an 
equivalent photon is incoherently launched into the cavity, either being radiated by the 
absorbing medium, or else, arriving from the vacuum beyond and being transmitted 
into the cavity. In any case, the probability that a photon with wavevector {Q, —q) 
and polarization e" arrives into the cavity from medium 2 with no phase relation to 
the lost photon is proportional to 1 — |r2|^. Similar statements apply to medium 1. 

From the previous discussion, it follows that in equilibrium the properties of the 
radiation field within the cavity V depend on the cavity walls only through their optical 
coefficients r"^. Equivalently, the cavity radiation is completely determined by the 
exact surface impedance Z^" defined through [fia x {fia x Ea)]'^ — Z^'^{na x HaY (sum 
implied), where fia denotes the outgoing unit normal of surface a, Ea and Ha are the 
total electric and magnetic fields at Za and fijV = x,y denote Cartesian coordinates 
along the walls. Thus, the electromagnetic radiation within the real cavity V would be 
identical to that within a fictitious cavity V' bounded by infinitely thin sheets at zi and 
Z2 , provided their refiection amplitudes r"^ are chosen to be equal to those of the walls 
of V. Their transmission amplitudes t"^ may then be chosen in order to guarantee 
energy conservation with no absorption whatsoever of electromagnetic energy (Fig. 
nib)). As there is no absorption in the fictitious system, there is no excitation of 
material degrees of freedom and the normal modes of the electromagnetic field form a 
complete orthogonal basis of the corresponding Hilbert space, allowing the use of well 
developed quantum- mechanical procedures for the calculation of the field properties. 
Contrariwise, in the real system the electromagnetic energy is absorbed, probably 
exciting electronic or vibrational transitions, so that the problem cannot be treated 
quantum mechanically without incorporating the electronic and/or vibrational degrees 
of freedom into the calculation, which would require in turn the use of a microscopic 
model of the material. 

In Fig. n^b) we have included two perfect mirrors at positions zq and z^ in order 
to quantize and count the normal modes of the system. The photons that are reflected 
at Zq and 23 and are transmitted back into V' mimic the photons injected into the real 
cavity V in order to restore thermal equilibrium, and in the limit i/, L777 — > 00 their 
phase is so large and so rapidly varying with lu that it effectively bears no relation 
with the phase of the photons lost from the cavity. 

Using the scattering approach we can treat dissipationless, homogeneous, 
isotropic, local, sharp media on the same footing as dissipative, inhomogeneous 
(layered systems, superlattices, photonic structures), chiral, spatially dispersive 
materials with a smooth selvedge. In particular, we can treat anisotropic systems. 
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3. Torque in ID 

We consider a finite beam propagating within V" along ±z, 

where the subindices r and / denote right and left moving contributions, Er and Ei 
are the corresponding amplitudes which we take as slowly varying functions of r, 
and Br and ei the polarizations within the x — y plane. To ensure that the field is 
divergence-less, an additional field contribution along z has to be added to (Q), 

A^(f, t) ^ (er ■ VEre'^"'-^*^ - • S/Eie"'^'' iz/q. (2) 



Expressions similar to Q and |2Jl may also be written for the magnetic field. The 
torque over medium 2 may be obtained by integrating the angular momentum flux 
M = r X T over a surface that surrounds it, where T is the electromagnetic stress 
tensor. Thus, 

r, = -i-Re J da[{fxE*)zEz + {fxB*)zBzl (3) 

where the integral is over the cross section of the beam. Notice that had we started 
our calculation with an infinitely extended plane wave, E^ and would have been 
zero, but the integral in Eq. (jS)) would have been over an infinitely extended surface, 
yielding an ill defined result. On the other hand, starting from Eq. Q we can take the 
limit of a plane wave, obtaining well defined expressions. The electromagnetic torque 
may be considered an edge effect that survives in the plane wave limit. Substituting 
and and similar expressions for the magnetic field i?, and after some manipulation, 
we obtain 

r. = --^Re(^ X dj*)z, (4) 

where A ^ oo is the cross sectional area of the wavefront. It can easily be checked 
that Eq. Q is consistent with the quantum mechanical view that each photon of 
energy Tiu) carries an angular momentum ±?i along ±z with speed c, according to its 
helicity. 

Now we consider one normal mode of the fictitious system with amplitude £q and 
frequency oj, E — £o(j){z)e~^'^'^ , where 

(f{z) = C^e"''' +D^e-''^\ {A = I,II, III) (5) 

is a spinorial normalized wavefunction with components 0^ {fi — x,y), C^, are 
distinct coefficients within each region A and q = lo/c. In the limit L/, L/// —^ oo, the 
electromagnetic energy U ^ [Li{\\C^\\^ + \\Dm^)+Liii{\\C'''W' + \\D''^\\^)]\£o\^A/8TT 
and the normalization condition 1 = [L/ ( 1 1 1 1 ^ + 1 1 1 1 ^ ) + L„ ^ ( j | C^^^ 1 1 2 + 1 1 1 1 2 )] 
are dominated by the large fictitious regions I and III, so we may identify U — 
A\£o\'^/8tt and solve for the amplitude \£o\'^ — STrft^huj/A in terms of the equilibrium 
photon occupation number /^^ — coth(/3?iti;/2)/2 at temperature fcsT = Thus, 
the contribution to the torque @ of one mode may be written as 
he 

Tz = --^L[(t>xdz(t>*y - (t)ydz(f>l + idz(l)y)(p* - {dz(Px)(l)*y]- (6) 

Now we label each mode by an index n and sum @ over n to obtain the total torque 

Tz ^ -he / dqfgc^6{q'^-ql)x[(j)nxdz(j)ly-(t)nydz(l3l^ + {dz(t)riy)(jjla:-i9z(IJnx)(IJn^^ 
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where we introduced the q integration and the Dirac's S in order to write the resuh 
in terms of the Green's function Gf^^{z,z') — J2,i'f'nfj,{z)4'nuW)/i'l^ +^V — 1n)i 



he f°° 
I" Jo 



(8) 



wheve A^^{z^ z') — [G'py(z, z') — Gy^(z', 2;)]/2i is the anti-Hermitian part of G^^. Here, 
?7 is a positive infinitesimal and we employed the identity 'k5(x) = —\Ta{x + iri)^^. The 
Green's function may be evaluated by solving {d1 + q^ + ir])Gfj,^{z, z') = S{z — z')5^^, 
subject to the appropriate boundary conditions. We write the solution in terms of 
the two homogeneous solutions u\{z) and v\{z) (A — 1,2) that satisfy the boundary 
conditions on the right and the left side of the system respectively, 

G{z,z')=u{z)[u'{z)-V{z')v-^{z')u{z')]-^e{z-z') 

- v(z)[v'(z') - u'iz')u-\z')viz')]-'eiz' - z), (9) 

where G(z,z') is a matrix with elements G^u{z,z'), u(z) and v(z) are matrices with 
columns ux{z) and v\{z) respectively, i.e., with matrix elements u^\{z) and v^\{z), 
u' and v' denote their derivatives with respect to their argument, and 9 denotes the 
Heaviside unit step function. 

For the case of uniaxial or orthorhombic slabs with normal-incidence reflection 
amplitudes r^^ and ry^ in the principal and Ha of the a-th slab we may write 

u(z) =R-Uo(z), v(z) = • vo(z), (10) 

where R and R^ are rotation matrices by angles 7/2 and —7/2 respectively and 



"0(^)^0 1 )^''''-'^{ rl je-^— ^ (11) 

1 A -^g(z-z,) .( fx, \ ,q(z-z,) 



vo(.)=^0 I je-^-^^ (12) 

are the solutions at the right and left sides of V' referred to the principal axes of the 
corresponding slab, which we assumed to be rotated by an angle 7 with respect to 
those of the opposite slab. Notice that R and R-^ act only on the first index /i of u 
and V. 

Substituting Eqs. ^-^^ in Eqs. ® and P we obtain a simple expression 

r°" ^ AriAra sin 27 6-2"-^ 

^"~~2^Jo AriAr2sin2 7e-2-i + (l-ri,r2,e-2«i)(l_r-i;,r-2j,e-2-i) ' 

where we also took the zero temperature limit fqc — 1/2 and we deformed the q 
integration path from the positive real axis toward the positive imaginary axis, as is 
usual. Here, k — q/i and we defined the anisotropy Ara ~ r^^ — Vy^ of the a-th slab. 
Eq. H13(l is the main result of the present paper. We have verified that it is equivalent 
to the result of Ref. jl3| . although in a much more compact form. 

4. Results 

As a first application of our result H13() we calculate the torque between two ideal 
perfectly reflecting mirrors covered by ideal perfectly absorbing polarizers, that is, 
we take Tx, = Tx^ = ±lj fyi = fy2 = 0. In this case, Eq. H13|l may be integrated 
analytically, yielding 

he 2 
= ——tan 7 log sin 7. (14) 
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Notice that the torque decays as 1/L, in analogy to the decay of the Casimir 

force between perfect mirrors in ID and the corresponding decay in 3D. In Fig. 
121 we show the torque as a function of the angle. As could have been expected, it is a 
periodic function of 7 with period tt. It is null when both polarizers are aligned, 7 = 0, 
corresponding to a stable equilibrium orientation, and when they are orthogonal to 
each other, 7 ~ ±tt/2, corresponding to unstable equilibrium. The inset shows that 
the slope of Tz{j) is singular at the stable equilibrium point. We remark that the 
torque is not simply proportional to sin 2j and therefore its extreme values are not 
at 7 = ±7r/4. However, it is odd-symmetric around 7 = 0. The maximum torque 
may be estimated as O.lhc/L; for example, at L = lOnm it is about lO^^^Nm. For 
dimensional reasons, in a full 3D calculation our result above would have to be scaled 
by A/L'^ multiplied by some dimensionless factor. 

If we replace the perfect mirrors above by lossy mirrors with reflection amplitude 
r, we can again obtain an analytical expression 

?ictan7 , I |2 2 \ /I r^ 

rz= ' log(l-|rpcos'7), (15) 

Fig. |31 shows that as \r\ diminishes Tz is reduced and becomes closer to a simple 
sinusoidal function « — fic\r\'^ sin 2^ /AttL. The inset shows that the singularity at 
7 = disappears when \r\ < 1. 

In Fig. ^we illustrate the torque between two identical dichroic mirrors relatively 
rotated by 7 = 7r/4 as a function of separation. Each mirror is characterized by a 
Lorentzian dielectric tensor with principal components 

(^) = 1 + ""f". , , (16) 

where is the frequency, r^^ the lifetime and w^i^p the intensity of a resonance 
of the a-th slab corresponding to polarizing field along the /Lt-th principal axis. In 
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-TT -37r/4 -7r/2 -7r/4 w/i 7r/2 Sn/A w 



angle a 

Figure 3. Torque between two lossy mirrors with reflection amplitudes |r| = 
0.6,0.7,0.8,0.9,1.0 covered by perfect polarizers as a function of the angle 7 
between corresponding principal directions. 
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this case the characteristic frequency LUp defines a characteristic lengthscale c/ujp. For 
separations much smaller than this distance the calculation is essentially non-retarded 
and the ID torque becomes constant, proportional to hujp. On the other hand, for 
larger separations we reach the retarded regime and the torque decays as hc/L, as for 
the ideal case. 

5. Conclusions 

Use of the scattering approach allowed us to obtain a simple expressions for the Casimir 
torque between anisotropic media in terms of their optical coefficients. Thus our results 
are applicable to arbitrary anisotropic materials and not only to semiinfinite, local, 
homogeneous ones. For instance, they may be readily applied to free standing or 
supported anisotropic films and to heterogeneous systems. We obtained analytical 
expressions for ideal systems which are the anisotropic counterparts to the ideal 
Casimir mirrors, and numerical results covering both the retarded and non- retarded 
regions for realistic dichroic systems with dispersive response functions. Our formalism 
has also permitted calculations of the torque between dissimilar materials, suggesting 
procedures to manipulate it, and it has been generalized to the full 3D case |18| . 
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